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Abstract 


A  new  measure  of  intrinsic  time  is  introduced  which  broadens 
the  endochronic  theory  and  lends  it  a  wider  predictive  scope.  Ideali¬ 
zed  plastic  models  are  shown  to  be  constitutive  subsets  of  the  general 
theory  and  the  phenomenon  of  yield  is  proved  to  be  a  consequence  of  a 
particular  definition  of  the  intrinsic  time  measure  in  terms  of  the 
plastic  strain  tensor. 

Various  versions  of  the  classical  plasticity  theory  are  shown 
to  be  asymptotic  cases  of  the  endochronic  theory.  In  particular  the 
kinematic  hardening  model,  the  isotropic  hardening  model  as  well  as 
their  combinations,  are  derivable  directly  from  the  general  theory.  In 
addition,  the  translation  vector  of  the  yield  surface  in  stress  space 
is  found  to  be  a  constitutive  property  given  by  a  linear  functional  of 
the  history  of  the  plastic  strain.  The  Prager  and  Ziegler  rules  are  im¬ 
mediately  obtainable  as  special  cases. 

It  is  also  shown  that  the  essential  features  of  plastic  re¬ 
sponse  under  conditions  of  "stress  reversals"  are  contained  in  the  con¬ 
stitutive  equation  whose  form  remains  invariant  of  the  deformation  his¬ 
tory.  It  is  believed  that  this  is  the  first  time  that  one  single  consti¬ 
tutive  equation  has  been  shown  to  predict  correctly  the  essential  features 
of  plastic  response  under  conditions  of  loading,  unloading  and  reloading. 
However,  when  the  intrinsic  time  measure  is  precisely  equal  to  the  norm 
of  the  increment  of  the  plastic  strain  tensor,  the  constitutive  unity, 
spoken  of  above,  no  longer  exists.  In  this  event  loading  and  unloading 
cases  must  be  distinguished. 


1.  Introduction 


The  concept  of  the  intrinsic  time  scale  was  introduced  as  a  proper  base 

of  measurement  of  the  memory  of  a  material  of  its  past  deformation  history, 

(1  2) 

leading  to  constitutive  theories  which  we  have  called  "endochronic" .  ’ 

The  case  of  strain  rate  independent  yet  history  dependent  materials  was 

dealt  with  at  length  in  previous  references  by  the  author  and  subsequently 

(  3“8  ) 

by  other  workers  in  the  field  whose  contributions  are  duly  referenced. 

Other  work,  in  other  directions,  involving  the  inelastic  behavior  of  metals 
is  currently  being  pursued  by  other  authors.  References  9,  10  and,  11  are 
typical  of  this  work.  In  the  papers  cited  heretofore,  we  introduced  two 
measures  of  intrinsic  time.  One  pertains  to  the  path  traced  by  the  deforma¬ 
tion  state  in  a  nine  dimensional  strain  space.  The  other  pertains  to  a 
stress  path  traced  by  the  state  of  stress  in  a  nine  dimensional  stress  space. 
If  we  denote  these  times  by  and  respectively,  then: 
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where  E  and  ■n  are  the  Green  strain  and  Piola  stress  tensors  respectively 
and  P  and  R  are  the  metrics  of  the  corresponding  spaces.  Specifically, 

L  |2t-  I*-  T..  (1.3) 
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where  T  is  the  Cauchy  stress  tensor.  Repeateci  indices  imply  summation  un¬ 
less  otherwise  stated. 

The  constitutive  equations  that  followed,  in  the  light  of  these  defini¬ 
tions,  were  obtained  from  thermodynamic  arguments.  In  particular,  the 
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internal  variable  theory  was  used  to  arrive  at  the  following  set  of  consti¬ 
tutive  equations,  which  pertain  to  isothermal  conditions, 


(1.4) 


(r  not  summed) 


(1.5) 


z  =  z(c) 


dz 

dC 


>  0 


(1.6) 


where  i|;  is  the  free  energy  density  and  are  the  internal  variables  of  the 
thermodynamic  system.  In  eq.  (1.5)  b^  is  a  dissipation  tensor  which  is  at 
least  positive  semidefinite  -  but  without  loss  of  generality  may  be  viewed 
as  positive  definite  since  otherwise,  if  =  0,  the  corresponding  r'th  Internal 
variable  may  be  eliminated  from  the  equations  as  a  consequence  of  the  re¬ 
sulting  condition: 


_  0  (1'7) 

^5r 


The  rate  of  irreversible  entropy  production  which  is  a  measure  of  the  in 
ternal  dissipation  is  obtained  from  the  equation 


(1.8) 


where  a  roof  over  a  symbol  indicates  differentiation  with  respect  to  z  and 
T  is  the  absolute  temperature.  Evidently  as  a  result  of  eqs,  (1.5)  and  (1,8) 
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Ty  =  q  •  b  •  q^ 

~  (r.  Slimmed) 

where  the  right-hand  side  of  eq.  (1.9)  is  an  inner  product,  clearly  a  posi- 
tive  scalar. 

Observation  of  accumulated  results  of  the  application  of  the  theory  to 
a  variety  of  histories  indicates  certain  broadly  consistent  trends.  The 
theory  is  evidently  simple,  versatile  and  has  powers  of  prediction  addition¬ 
al  to  the  theories  of  plasticity  of  the  classical  type  provided  there  are  no 
reversals  in  the  rate  of  stress.  In  particular,  in  one  dimension,  the  slope 
of  unloading  at  a  point  as  the  uniaxial  stress  strain  curve  was  predicted  by 
the  "linear"  version  of  the  theory  to  be  2E^  -  where  E^  is  the  elast* 
modulus  and  the  tangent  modulus.  This  Is  an  over-estimate  of  the  normaUy 

observed  unloading  slope  which  is  close  to  E^. 

The  purpose  of  this  paper  is  to  eliminate  this  deficiency  of  the  endo- 

chronic  theory  by  Introducing  a  measure  of  intrinsic  time  which  is  more 
closely  representative  of  the  dissipation  properties  of  metals.  The  under¬ 
lying  cause  of  this  discrepancy  is  of  a  thermodynamic  nature  and  has  to  do 
with  the  fact  that  the  theory,  heretofore,  predicts  a  rate  of  dissipation 
during  loading  which  is  identical  to  that  which  occurs  at  the  onsent  of  un¬ 
loading  . 

In  this  paper  we  introduce  a  new  concept  of  intrinsic  time  ?  which  cor¬ 
rects  both  these  deficiencies.  In  one  (axial)  dimension  we  stipulate  that 

dC=  Ids-kf  1  (1-10 

o 

where  k  is  a  positive  scalar  such  that  0  <  k  <  1  and  E^  is  the  elastic  modu¬ 
lus.  If,  at  the  onset  of  unloading,  we  denote  the  unloading  slope  by  E_ 
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and  the  rate  of  dissipation  by  y_,  then  it  is  shown  that  E_  tends  to  and 
Y  tends  to  zero  as  k  tends  to  vinity. 


In  three  dimensions  we  form  a  strain-like  tensor  6^^  given  by  eq.  (1.11) 


9.  . 

1] 


(1.11) 


where  ({>  is  a  positive  definite  symmetric  fourth  order  material  tensor.  Deal¬ 
ing  strictly  with  isotropic  materials  where  (j)  is  isotropic  we  proceed  to  de¬ 
fine  a  deviatoric  strain-like  tensor  by  the  equation: 


e  •  •  ~  7r~  s .  . 
i:  2y  1] 


(1.12) 


where  e  and  s..  are  the  deviatoric  strain  and  stress  tensors  respectively. 
13  13 

We  also  define  a  hydrostatic  strain-like  tensor  6^^  by  the  equation 
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The  relation  between  k  and  k  on  one  hand  and  the  components  of  (()  on  the  other 

o  1  ; 

is  shown  in  Section  6.  Evidently,  2^^  and  9^^  are  the  deviatoric  and  hydro- 

Static  components,  respectively,  of  ^^j* 

A  hydrostatic  intrinsic  time  measure  and  a  deviatoric  counterpart 

are  now  defined  by  the  eq.*s 


"4  = 

"  -01  "  2l; 

.  d  2- 

1  ^13 

(1.14a) 

^4- 

^10 

"-11 

,  d2..  • 

!  13 

(1.14b) 

where  k  are  material  parameters  which  provide  for  a  coupling  between  hydrO' 
rs 

static  and  shear  response. 

The  resulting  constitutive  equation  for  isotropic  materials  is  given 
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by  eq. ’s  (6.8)  and  (6.9).  These  are  identical  in  form  to  those  of  the  simple 
endochronic  theory  (where  both  and  k^  are  equal  to  zero)  but  predict  be¬ 

havior  that  is  far  closer  to  that  of  metals,  when  k^  and  k^  are  close  to 

unity. 

The  very  significant  particular  case:  k^  =  1,  z^  =  0. 

This  is  the  case  of  elastic  hydrostatic  response  and  the  intrinsic  time 
measure  dC  is  equal  to  the  norm  of  the  plastic  strain  tensor.  The  resulting 
constitutive  equation  is  eq.  (6.36)  or,  equivalently,  (6.38).  With  the  aid  of 
these  equations  the  following  propositions  are  proved  in  Section  6  of  the  text. 

(i)  A  spherical  yield  surface  in  deviatoric  stress  space  exists. 

(ii)  If  =  0  in  eq.  (6.36)  and  f(C)  increases  monotonically  with 
then  the  classical  theory  of  plasticity  with  isotropic  hardening  follows.  The 
increment  of  plastic  strain  is  shown  to  be  normal  to  the  yield  surface. 

(lii)  If  #  0  (in  which  case  p^^  >  0)  and  f(?)  =  1,  then  the  spheri¬ 
cal  yield  surface  translates  in  deviatoric  stress  space  and  a  theory  of  kine¬ 
matic  hardening  results  with  a  general  rule  which  contains  Prager  s  and 
Ziegler’s  rules  as  special  cases. 

(iv)  If  p^  5^  0  and  f(  C)  is  a  monotonically  Increasing  function  then 
the  yield  surface  translates  and  expands  simultaneously  according  to  rules 
inherent  in  the  theory. 

(v)  The  above  are  true  in  the  case  of  softening,  instead  of  hardening, 
when  f( O  is  a  monotonically  decreasing  function. 

Remark:  No  yield  surface  exists  if  <  1.  Thus  the  classical  theory 


/ 
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of  plasticity  is  the  "boundary"  of  the  endochronic  theory. 


2.  The  Rate  of  Entropy  Production 


For  the  purposes  of  concentrating  on  the  physical  aspects  of  the  prob¬ 
lem  and  to  rid  our  argument  of  superfluous  analytical  complications,  we 
limit  the  discussion  in  this  Section  to  a  one  dimensional  stress  field  and 
to  theories  with  one  scalar  internal  variable.  To  this  end  we  let 


=  3sE  (e  -  q)" 


(2.1) 


a  =  ||-  E(e  -  q) 


(2.2) 


or. 


If 


a  =  bq 


=  0 


(2.3) 


(2.4) 


where  e  is  the  uniaxial  strain  and  a  the  appropriate  axial  stress.  No  limi¬ 
tation  to  small  deformation  is  therefore  implied  except  in  so  far  as  the 
validity  of  the  algebraic  forms  of  eqs.  (2.1)  and  (2.3)  are  concerned.  It 
follows  from  eq.  (1.9)  that 


dz  b 


(2.5) 


A  typical  unloading  response,  given  by  the  linear  version  of  the  endo¬ 
chronic  theory,  using  eqs.  (2.2)  and  (2.3)  is  shown  in  Figure  1.  The  slope 


at  the  point  of  unloading  (point  A)  is  2E 


E^,  in  variance  with  the  observed 
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slope  which  is  in  the  vicinity  of  E^.  The  reason  for  the  discrepancy  may 
well  lie  in  the  expression  for  the  dissipation,  i.e.  eq.  (2.5).  The  in¬ 
ference  is  that  the  rate  of  dissipation  at  A  is  the  same  whether  one  pro- 
ceeds  to  load  or  unload  at  A. 

This,  unfortunately,  is  not  in  accord  with  experimentally  established 

mechanisms  of  plastic  flow.  For  instance,  in  single  crystals  the  latter 

occurs  during  loading  by  cumulative  slip  along  planes  between  which  the 

(12  ) 

crystal  suffers  only  elastic  distortion.  During  initial  unloading, 

back  slip  is  minimal  and  the  crystals  change  shape  mainly  (if  not  totally) 
by  elastic  recovery  of  the  parts  of  the  crystal  between  slip  planes .  This 
indicates  that  the  rate  of  dissipation,  if  not  zero  during  unloading, 
should  at  least  be  much  less  than  its  counterpart  during  loading. 

The  endochronic  theory  cannot  predict  such  a  result  unless  the  incre¬ 
ment  of  intrinsic  time  corresponding  to  a  decrement  de_  is  much  less  than 
that  due  to  an  increment  de_^  such  that 

=  \^^+\ 

If  indeed 

dC  =  k|de  1 

where 

k  =  k(e,  a,  C) 

a  prediction  of  this  type  of  dissipation  behavior  is  impossible.  This 


(2.6) 


(2.7) 
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observation  lies  at  the  root  of  the  difficulty  that  one  encounters  in  using 
the  endochronic  theory  to  describe  unloading. 

Interestingly  enough,  quasilinear  versions  of  the  theory  have  been 
proposed^ which  correct  the  slope  of  the  unloading  response.  However, 
the  thermodynamic  relation 

dy^  =  dY_  (2.9) 

* 

in  obvious  notation,  remains  intact  and  equally  inapplicable  to  metals. 


*To  avoid  this  difficulty  other  authors^  ’  ^  have  proposed  one  type  of 

constitutive  equation  for  loading  and  another  for  unloading.  Here  we 
show  that  the  difficulty  is  circiamvented  without  this  dichotomy. 
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The  situation  is  similar  when  many  internal  variables  are  used.  In 
this  case  the  expression  for  p  is  a  generalization  of  eq.  (2.1): 


♦  =  I  I  "^r  ■  '*r'‘ 

r 


(2.10) 


correspondingly , 


o  =  I  (e  -  q^) 

r 


(2.11) 


+  b  =  0 


9qr 


(2.12) 


T  ^  =  -y  i  . 


(2.13) 


This 


is  last  equation  in  conjunction  with  eq.  (2.12)  shows  that 


T  ^  =  y  1.  (^) 
^  dz  ^  b. 


(2.14) 

r  -r  'V 

The  important  conclusion  from  the  right-hand  side  of  eq.  (2.14)  is  that  ^ 
is  a  function  of  state,  i.e.,  of  the  state  variables  e  and  q^.  Thus  if 


dz^  =  dz_ 


(2.15) 


when  lde|^  =  ldel_.  then 


dY^.  =  dY_ 


(2.16) 


as  before. 
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From  the  above  discussion  it  is  fairly  conclusive  that  the  unloading 
behavior  of  metals  will  be  predicted  more  precisely  if  dz  <  dz^  at  the 
point  of  unloading,  for  the  same  value  of  |de|. 


3.  Physical  Considerations 

The  essential  premise  of  the  endochronic  theory  is  that  the  intrinsic 
time  increment  is  a  measure  of  irreversibility.  For  instance,  reference  to 
eq.  (2.3)  shows  that 

Thus  if  in  the  course  of  a  change  in  the  thermodynamic  state  of  a  system  dz 

is  zero,  then  it  follows  from  eq.  (3.1)  that  dq  is  also  zero,  in  which  case 

(13  ) 

the  system  is  behaving  reversibly.  Since  at  the  point  of  unloading,  the 

process  of  unloading  borders  on  the  reversible,  though  it  may  not  be  exactly 

dz 

so,  one  is  led  to  conclude  that  in  one  dimension  the  ratio  i'--'  T  -  must  be 

W\ 

close  to  zero.  Furthermore,  since  during  elastic  unloading  in  one  dimension 
the  strain  decrement  is  governed  essentially  by  the  relation 


(3.2) 


one  is  led  to  the  conclusion  that  a  definition  of  intrinsic  time  which  is 
more  closely  representative  of  unloading  behavior  should  be  given  by  the  re- 
lat ion 


d;  =  |de  -  k  1^1 
0 


(3,3) 
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where  k  is  a  scalar  which  in  the  case  of  metals  is  close  to  and  possibly 
equal  to  unity  and  certainly  non-negative .  Thus : 

0<k<l.  (3.4) 

Note  that  except  for  a  scalar  multiplier  this  definition  of  intrinsic  time 
reduces  to  our  previous  one ,  ^  ^  when  k  =  0 . 

It  will  be  shown  presently  that  as  k  approaches  unity  the  constitutive 
behavior  as  depicted  by  a  single  internal  variable »  approaches  that  of  an 
elastic  perfectly  plastic  solid. 

4.  Application  of  the  Theory  to  One  Dimension 

We  proceed  to  apply  the  theory  to  those  materials  which  do  not  harden 
cyclically.  In  this  case  z  =  ?;  • 

Elimination  of  the  internal  variable  q  from  eqs.  (2.2)  and  (2.3)  leads 
to  the  following  constitutive  equation  in  one  dimension: 

[de-k^l  a=b(de-  )  .  (4.1) 

o  o 

We  shall  proceed  to  study  materials  for  which  de  >  ^  .  In  this  case  the 

o 

sign  of  de  -  k  . is  the  same  as  the  sign  of  de. 

ilj 

o 

We  shall  consider  the  cycle  where  loading  begins  at  (  a, e)  =  (0,0)  and 
continues  into  the  plastic  range  where  unloading  begins,  followed  by  com¬ 
pression  into  the  plastic  range  and  unloading  to  zero  stress. 

In  Figure  3,  we  show  a  plot  of  ^  vs  /  [  de]  for  a  perfectly  plastic 
solid  for  the  purposes  of  comparison  with  the  predictions  of  the  constitu¬ 
tive  equation  (4.1). 
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In  the  course  of  the  paths  OYA  and  CDEA  of  the  ideal  solid  of  Figure 
3,  eq.  (4.1)  will  read 

d.  -k|2^=b(dc-|i)  .  (“.2) 

o  0 

It  follows  from  this  equation  that 

L_  ^  ^  ~  (4.3) 

E  de  ~  1  -  k(a/b) 

0 

As  a  increases  the  numerator  will  vanish  before  the  denominator  at  b  =  o 
since  k  <  1.*  Furthermore,  it  can  be  easily  verified  that  as  k  -*■  1,  the 
slope  4^  of  the  stress  response  approaches  that  of  a  perfectly  plastic 
solid  given  in  Figures  2  and  3,  for  the  paths  OYA  and  CD. 

In  the  course  of  the  unloading-compression  path  ABC  of  Figure  2 ,  eq. 
(4.1)  will  read 

1_  ^  -  1  t  (g/b)  (4.4) 

E  de  “  1  +  k(o/b) 

0 

which  predicts  a  behavior  which  will  also  tend  to  that  of  the  perfectly 
plastic  solid  as  k  1. 

Also  of  interest  is  that  eqs.  (4.3)  and  (4.4)  give  stress  responses 
in  tension  and  compression  which  are  mirror  images  of  each  other  as  one 
can  readily  verify  by  substituting  minus  c  for  a  in  eq.  (4,4), 

Integration  of  the  constitutive  equation  (4.1)  for  the  two  conditions 
de  >  0  and  de  <  0  gives  the  following  two  equations  respectively  for  the 

*  The  case  of  k  =  1  must  be  treated  as  a  limiting  case  otherwise  one  is 
lead  to  erroneous  results. 
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uniaxial  stress  response; 

de  >  0; 


TT  = 

b  o 


0-0 

k  -7-^  -  (1-k)  log 
b 


b-g 

b-o 


(4.5) 


de  <  0: 


^0  (e-E  )  = 
b 


k 


0-0 

+  (1-k)  log 
b 


b+0 

b+o 


(4.6) 


where  (e  ,o  )  is  the  initial  point  on  the  path  of  the  appropriate  solution. 

0  0 

The  essential  constants  of  eq.  (4.5)  are  the  initial  modulus  and  the 
ultimate  stress  b.  In  Figure  4  we  give  a  stress  response  curve,  typical  of 
pure  aluminum  where  =  9.1  x  10®  Ib/in^  and  b  =  16.24  x  10®  Ib/in^  It 
is  quite  apparent  that  the  material  response  tends  to  that  of  an  elastic  per¬ 
fectly  plastic  solid  as  k  ^  1. 

This  result  is  quite  obviously  a  major  development  in  the  evolution 
of  the  endochronic  theory. 


5.  The  Question  of  Dissipation 

The  difficulty  that  one  encounters  in  calculating  the  dissipation, 
using  the  older  version  of  the  endochronic  time  measure  given  by  eq.  (1.1) 
has  also  been  circumvented.  The  following  expressions  are  easily  obtained 

for  dy^  and  dy  : 

Tdv  =  (A -iV  1-^4 
*  l-k(r) 
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dc 


(5.2) 


If  one  takes  limits  of  dy  in  their  proper  sequence,  one  finds  that 


Lim  (Lim^^j^  dy^)  =  b|dc| 
k->l 


(5.3) 


Lim  (Lim  dy  )  =  0 

,  .  o->-D 

k->l 


(5.4) 


Eqs,  (5.3)  and  (5.4)  give  the  precise  dissipation  properties  of  an  elastic 
perfectly  plastic  solid  in  the  plastic  state.  In  Figure  5  ^  is  plotted 

versus  stress  in  the  case  of  a  solid  where  k  =  .99  on  a  semilogarithmic 
scale.  Note  that  as  the  deformation  increases  the  divergence  between 
and  increases  the  limiting  ratio  being  equal  to  2  x  10^. 


6.  Generalization  to  Three  Dimensions  and  n  Internal  Variables 

We  define  the  strain-like  tensor  0^^  by  the  following  equation: 


®ij  "  ^ij  ~  '*‘ijKl!. 


(6.1) 


where  ^  is  a  positive  definite  symmetric  fourth-order  tensor.  In  the  case 
of  isotropic  materials,  $  is  of  the  form: 


d  =  5.,  6  „  d>  +  (6.  6.„  +  d.„6  .)  . 

^ij<«.  ij  kJI  IK  Kj 


(6.2) 
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In  this  event,  0  may  be  decomposed  into  its  deviatoric  and  hydrostatic  parts 
as  follows: 


KIC  ij 


(6.3) 


where 


k 

o 

®kk  ”  ^kk  ”  3K^  kk 


(6.4) 


(6.5) 


In  eq.  (6.5)  e  and  s  are  the  deviatoric  parts  of  e  and  a  respectively,  jj^ 

and  K  are  the  elastic  shear  and  bulk  moduli  of  the  material  respectively 
o 

and  the  constants  k^  and  k^  are  related  to  and  by  the  relations 


2^1 


k 

o 


various  degrees  of  generality  are  now  possible.  For  instance,  mindful 

of  the  consequence  that  in  the  case  of  isotropic  materials,  undergoing  small 

(2) 

deformation,  the  deviatoric  and  hydrostatic  responses  are  separable,  we  may 
define  a  hydrostatic  intrinsic  time  measure  dc^^,  where 


00  kk  ll  01  13  13 


(6.6) 


and  a  deviatoric  intrinsic  time  measure  dc^  where 

<^'■1  =  «10  ‘'\k  ““ti  *  “U  • 


(6.7) 
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The  hydrostatic  measure  would  then  be  the  appropriate  one  to  use  in  the 
hydrostatic  constitutive  response,  while  dc^  would  be  used  in  the  deviatorxc 
constitutive  response.  In  this  case  (i.j  =  1,2)  is  a  matrix  of  nondimen 

sional  scalars. 

Using  the  formulation  of  Reference  1  as  a  point  of  departure,  we  may 
then  write  the  rate  independent  (plastic)  response  of  metals,  in  the  small 
deformation  region  as  follows; 


§  =  2 


3e 


dzj 


(6.8) 


and 


a 


KK 


dz  ’ 
H 


(6.9) 


where 


and 

- 

« ■ 


The  functions  and  f  are  both  non-negative . 
D  n 


(S'. 10) 


(6.11) 
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6.1  Particular  Cases 

We  wish  to  examine  in  greater  detail  the  case  where  1  and, 

K  and  K  „  are  equal  to  zero.  We  also  set  dC„  identically  zero  there- 
ol  10  ^ 

by  dealing  with  a  material  with  an  elastic  hydrostatic  response  and  no 
coupling  between  deviatoric  and  hydrostatic  behavior.  Now  the  constitu¬ 
tive  description  of  the  material  simplifies  considerably  so  that 


a 

KK 


3  K  e 

KK 


(6.12) 


z  3e 

s  =  2  f  ii(z  -  z') 
o 


(6.13) 


where  we  have  temporarily  dropped  the  suffix  D  in  since  no  confusion 
is  likely  to  result.  Also  we  have  set 


=  dOj. 


(6.14) 


®ij 


-  k 


Ji 

1  2y^ 


(6,15) 


and 


dz  =  ,  f(C)  >  0  . 

The  constant  in  eq.  (6.7)  does  not  appear  in  eq.  (6.14)  since  it  is  now 
an  immaterial  constant  and  can  be  appended  to  the  function  f(c).  Such  re¬ 
dundancies  in  the  number  of  constants  that  describe  the  constitutive 
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properties  of  a  metal  are  present  in  the  theory  and  can  be  dealt  with  as 
they  arise  without  the  need  for  further  comment. 

6.1.2  Uncoupled  Hydrostatic  and  Deviatoric  Responses 
The  nature  of  the  matrix  merits  further  discussion.  For  in¬ 
stance,  if  it  is  diagonal  then  there  is  no  coupling  between  the  deviatoric 
and  hydrostatic  response.  In  this  case  a  hydrostatic  strain  history  will 
have  no  effect  on  the  shear  behavior  and  vice  versa.  On  the  other  hand 
the  theory  does  admit  such  an  effect  if  is  not  diagonal.  At  this  time 
we  do  not  have  enough  evidence  to  comment  on  the  symmetry,  or  asymmetry, 

of  K.j, 

Of  particular  importance  are  the  cases  where  the  deviatoric  and  hydro¬ 
static  responses  are  uncoupled  and  either  k^,  k^  or  both  are  equal  to  unity. 
We  wish  to  discuss,  at  length,  the  case  of  the  shear  response;  the  arguments 
will  apply  equally  well  to  its  hydrostatic  counterpart.  With  reference  to 
eqt  (6,8)  one  obtains  the  following  relation 


s 


da 


(6,17) 


where 


y(z)  G(z)  ;  G(0)  =  1  , 


C6.18) 


and  p(z)  is  related  to  G(z)  by  the  integral  equation 


p(z)-k^f  p(z.2')  g,  dz»  =  G(z)  . 


(6.19) 


The 


solution  of  this  equation  in  the  case  of  two  internal  variables  is  given 


in  Appendix  A. 
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More  simply,  the  following  algebraic  relation  exists  between  the  Laplace  trans¬ 
forms  p  and  G  of  p  and  G,  where  p  is  the  Laplace  transform  variable: 

G 


P  = 


(6.20) 


1-k^pG 


This  specific  form  of  the  theory  has  the  characteristic  feature  that  it  leads 

to  a  constitutive  equation  which  gives  the  deviatoric  stress  response  to  the 

history  of  the  deviatoric  strain- like  tensor  Z  in  terms  of  the  path  in  Z 

space.  We  point  out  however  that  there  is  no  specific  connection  with  an 

a  priori  existence  of  a  yield  surface.  Furthermore  there  is  no  dichotomy 

in  the  constitutive  representation  of  the  loading  and  unloading  responses. 

Equations  (4.5)  and  (4.6)  are  derivatives  of  one  and  the  same  eq.  (4.1). 

It  is  of  interest  to  note  that  d2  is  akin  to  and  is  equal  to  the  de- 

viatoric  "plastic  strain"  increment,  in  classical  plasticity  terms,  if  k^  is 

ds 

unity  and  the  "elastic  strain"  increment  is  defined  as  ^  .  Moreover,  we 
will  proceed  to  show  that  if  k^  =  1,  then  a  yield  surface  exists.  The  proof 

will  apply  strictly  to  the  case  where  k^  =  1. 

To  this  end  we  invoke  the  result  of  Reference  1  according  to  which: 


G(z) 


I 


r=l 


G 

r 


'KX  Z 

r 


e 


C6,21) 


where  G  are  positive j  also,  are  all  positive  with  the  exception  of  a 
which  may  be  zero.  Furthermore,  since  G(0)  =  1  it  follows  that 


n 


I 

r=l 


G 

r 


1 


C6.22) 


Let  G(p)  be  the  Laplace  transform  of  G(z),  Then  as  a  result  of  eq,  (6,21) 
G(p)  is  of  the  form  _ _ 


*For  explanatory  notes  on 


this  treatment  see  Appendix  C. 


20 


=  Iff} 


(6.23) 


where 


P  =  G^(p  +  a^)  ...  (p  +  a^)  +  (p  +  cij)  ...  (?  +  «*„) 

+  (p  +  a^)  ...  (p  +  «n.i)  Gjj 


(6.24) 


Q(p)  =  (p  +  a^)  (p  +  “2^  *"  “n^  ' 


(6.24) 


Evidently,  P  and  Q  are  polynomials  of  order  (n-1)  and  n  respectively,  Further¬ 
more,  the  coefficient  of  the  leading  term  of  P  is  equal  to  unity.  It  follows 
from  eqs,  (6.20)  and  (6.23)  that 


“‘P’  =  R(p) 


(6.25) 


where 


R(p)  =  Q(p)  -  p  P(p)  . 


(6,26) 


It  may  be  shown  by  direct  computation  that  R(p)  is  a  polynomial  of  order 

n-1  " 

n-1 ,  such  that  the  coefficient  of  p”  is  exactly  the  sum  ^  a  G  .  Hence 

P(  ^  ^ 

the  ratio  is  not  function-like,  in  the  sense  of  Mikusinski,  in  that 

R(p) 

it  contains  a  delta  function  of  strength  - - -  .  In  fact  if  -6^  are  the 

y  a  G 
r=i  r  r 


zeros  of  R,  then  p(p)  may  be  written  as 
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n-1  R 


p(p)  = 


+  I 


U  0  ^ 

r  r 

r=l 


(6.27) 


since  p(p)  - 


n 

y  a  G 

u  r» 


is  the  ratio  of  two  polynomials,  the  numerator  being 


r=l 


r  r 


of  degree  n-2,  and  the  denominator  R(p)  of  degree  n-l.  It  is  shown  below  that 
if  the  absolute  values  of  the  zeros  of  Q(p)  are  ordered  in  the  sense  that 


<  02  <  Og  ...  <  o^ 


then  6  must  always  satisfy  the  inequalities; 
r 


“l 


o-  <  <  o. 


e  <  o 
n-l  n 


and  are  therefore  all  positive.  The  proof  is  elementary.  Evidently 
R(-a^)  =  ***  ° 

R(-a2)  =  “2^  ^“n~  ^  °  ' 


Similarly, 

RC-Og)  <  0 

and  generally  R(p)  alternates  in  sign  at  the  zeros  of  Q(p).  It  follows  that 

R(p)  must  vanish  at  the  points  p  =  -6^,  P  =  -82  •••  P  =  ”^n-l 

all  positive  and  are  bounded  from  above  and  below  by  o^  according  to  the 
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inequalities  given  above.  Finally,  since  R(p)  is  a  polynomial  of  degree 
n-1  it  can  have  no  other  zeros.  This  completes  the  proof. 

The  values  of  the  residues  are  found  from  the  formula 


R  = 
r 


dR 


(6.28) 


R 


where  R  (p)  def  ^  .  For  a  detailed  calculation  see  Appendix  B  where  it  is 
-dp 

shown  that  the  residues  R^  are  all  positive. 

Clearly,  as  a  result  of  eq.  (6.27) 


6(z) 


n-1  -6  z 

r 


p(z)  =  ♦  I  R.  e 


(6.29) 


n 


y  o  6 
r  r 

r=l 


r=l 


which  we  write  as 


p(z)  =  p  6(z)+  p^(z) 


(6.30) 


where  p-^(?;)  is  composed  of  a  finite  sum  of  exponential  terras, 
Evidently,  it  follows  that 


d(i  -z  3<Z 


(6  .31) 


Now,  as  a  consequence  of  eq.  (6.31),  at  z  =  0: 


s  =  2vi 

-  0  0 


da 

dz 


z=0 


(6  .  32) 


Also  from  eq.  (6.15)  (and  for  =  l),.the  condition  a  =  0  (and,  therefore, 
z  =  O)  gives  the  relation 
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s'=  2jj  e 
0 


(6.33) 


Equation  (6.33)  merely  attests  to  the  fact  that  while  z  =  0  the  deformation 
process  is  reversible  and  therefore  the  deviatoric  stress  response  is  elastic. 
See  Sec. 7.  It  is  also  of  interest  that  at  z  =  0,  ^  is  indeterminate  and 
can  take  any  value  consistent  with  eq,  (6.33).  Specifically,  eqs.  (6.32) 
and  (6.33)  combine  to  give 


(6.34) 


da 

However,  at  z  >  0  the  derivative  exists  and  specifically  at  z  =  0+  one 

da 

obtains  the  limit  of  by  approaching  z  =  0  from  the  right.  In  fact, 

QZ 

the  point  0+  is  the  point  of  deviation  from  elastic  response  or  the  yield 
point.  At  this  point 


in  accordance  with  eq.  (6.14).  Therefore,  applying  eq.  (6.33) 

In  conclusion^  in  the  process  of  monotonic  loading,  while  |s|  ^  z  =  0 
and  eq*  (6.33)  applies  and  the  material  response  is  elastic.  However, 
when  [s]  >  s^  9  then  eq.  (6.31)  applies  and  the  response  is  no  longer  elas¬ 
tic.  Of  course,  s^  is  the  YIELD  stress  and  eq.  (6.35)  is  the  VON  MISES 
yield  criterion. 
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In  particular,  as  a  result  of  eq.  (6.35),  the  constitutive  relation 
(6.31)  may  be  vfritten  more  succinctly  as 


s 


(6.36) 


where  s°  =  2vi  p  and  has  the  physical  significance  of  an  initial  yield  stress. 
Y  0  0 

We  shall  show  presently  that  eq.  (6.36)  yields  very  rich  results  and  re¬ 
veals  important  characteristics  of  plastic  behavior  that  first  appeared  as 
assumptions  or  conjectures  in  the  classical  theory  of  plasticity. 

To  this  end,  let  the  Integral  on  the  right-hand  side  of  eq.  (6.36)  be 


denoted  by  r,  i.e.,  set: 


r  =  2v^  /  p^(Zp  -  zy  ^  dz'^ 


(6.37) 


Equation  (6.36)  then  reads, 


dj 


(6.38) 


We  recall  that  =  1  and  therefore  d^  is  exactly  equal  to  the  increment  of 

plastic  strain.  We  wish  to  examine  eq.  (6.38)  in  the  specific  case  where 

and  C  are  uncoupled, 

H 


d^^  =  d2  *  d^ 
D  ~  ~ 


(6.39) 


and  the  hydrostatic  response  is  elastic,  i.e.,  dz^^  =  0,  in  which  event  eq. 
(6.9)  gives  rise  to  eq.  (6.12),  i.e. 


a  =  3Ke 

KK  KK 


in  our  previous  notation. 


(6.40) 


/ 
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Case  (i):  f(i;)  =  1 

Equations  (6.38)  and  (6.39)  give  rise  immediately  to  the  following  two 
results : 


(6.41) 


d^  =  —  (s  -  r)  de^  (6.42) 

It  is  clearly  obvious  that  if  r  =  0  then  eq.  (6.42)  is  that  of  an  elas¬ 
tic  perfectly  plastic  material  with  a  Von  Mises  yield  criterion.  In  this 
case 


(6 .42a) 


On  the  other  hand,  if  r  i  0,  then  eq.  (6.41)  shows  readily  that  case  (i) 
corresponds  to  kinematic  hardening.  This  relation  is  in  fact  the  equation  of 
a  hypershere  in  deviatoric  stress  space.  It  also  represents  the  equation  of 
a  circle  in  principal  stress  space.  In  either  case  s°  is  the  radius  of  the 
hypershere  (circle)  and  r  is  the  radius  vector  which  connects  the  origin  of 
the  stress  space  to  the  centre  of  the  hypersphere  (circle) .  The  yield  surface 
is,  therefore,  a  translating  spherical  (circular)  surface.  Eq.  (6.42)  shows 
that  the  increment  in  plastic  strain  is  normal  to  the  yield  surface.  These 
results  are  shown  diagrammatlcally  in  Figure  S .  Note  however  that  whereas  in 
the  classical  theory  of  plasticity  the  concept  of  kinematic  hardening  was  a 
conjecture,  here  it  is  a  derived  result > 

Furthermore  5  we  feel  it  Important  to  emphasize  that  in  the  classical 
not  known  how  the  surface  translates,  i.e.,  it  is  not  known 


theory  it  is 


a  priori  how  r  depends  on  the  history  of  loading  or  plastic  strain 
stance,  Prager  assumed  that 
dr  =  cd  ^ 

where  c  is  a  material  constant.  Ziegler  suggested  that 


For  in- 


(6.43) 


where  dy  is  a  positive  quantity  not  specified.  We  point  out  that  eq.  (6,43) 
is  equivalent  to  eq.  (6.44)  if  dy  is  proportional  to  d^jj  as  eq.  (6.42)  readily 
Indicates.  However,  Prager’s  (or  Ziegler’s)  rule  is  a  particular  case  of  the 
present  theory  as  is  pointed  out  in  the  following  Remark. 

Remark  1:  Prager 's  rule  of  kinematic  hardening  is  a  particular  case  of  the 
present  theory  obtained  from  eq.  (6.37)  but  setting 


Pl(^p)  =  constant 
in  which  event. 


(6.45) 


(6.46) 


This,  as  is  well  known,  is  called  linear  hardening  in  the  sense  that  the 
stress  strain  curve  in  shear  is  linear  beyond  the  onset  of  yield. 

In  this  same  vein,  if  P^CC^)  is  not  a  constant  but  consists  of  a  single 

exponential ,  i . e . , 

Pl(Cp)  =  P^e"“^D  (6.47) 

then  it  is  easily  shown  as  a  result  of  eq.' s  (6.37)  and  (6.38)  that 


(Jp  =  (J2;(s  -  Br) 

„  o  „  ~ 


where 


as. 


6  =  1  + 


2H>Pi 


(6.48) 


(6.49) 
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This  is  a  new  type  of  kinematic  hardening.  Note  that  the  yield  surface  no 
longer  translates  along  the  outward  normal  at  the  extremity  of  the  stress 
vector  but  in  a  direction  which  is  skew.  The  skewness  depends  on  the  value 
of  6. 

Of  course  the  above  are  particular  cases.  In  general  the  translation 
vector  is  determined  by  eq.  (6.37),  in  terms  of  a  convolution  product  which 
involves  the  material  function  pj^(Sp).  Within  the  assumptions  of  the  present 
case,  it  is  important  to  observe  that  can  be  determined  by  a  simple 

shear  experiment.  In  fact  one  can  show  readily  from  eq.  (6.37)  that 


dx 


.T  > 


(6.50) 


where  x  is  the  shear  stress  and  the  tensorial  plastic  shear  strain  compo¬ 
nent  . 

Remark  2:  The  function  is  determinate  from  a  single  monotonic  shear 

test . 

Remar^k  3:  The  mode  of  translation  of  the  yield  surface  is  determined  from 
a  single  monotonic  shear  e5q)eriment. 

Case  (ii):  f ( g)  monotonically  increasing 

The  counterparts  of  eq.*s  (6.41)  and  (6.42)  are  now  the  following: 

1  Is  -  r] I  =  f( ?q)Sy  (6.51) 


de 


(s  -  r)  dCp 


(6.52) 


*for  all  histories 
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Clearly,  if  r  =  0,  then  eq.  (6.52)  is  that  of  a  plastic  material  with 
isotropic  hardening  and  a  Von  Mises  yield  criterion.  In  this  case 


(6.52a) 


Sy 


On  the  other  hand  if  r  ^  0  then  eq.  (6.52)  shows  that  the  yield  surface  now 
expands  as  well  as  translates.  The  increment  of  plastic  strain  is  still  nor¬ 
mal  for  the  yield  surface.  The  translation  vector  r  is  still  given  by  eq. 
(6.37).  We  shall  not  go  into  the  details  of  determining  f(C)  and  P3^(^d)  ^ 
this  case,  at  this  juncture.  The  reader  is  referred  to  Appendix  D. 

The  constitutive  equation  for  a  >  Oy  for  the  hydrostatic  response  is 
identical  in  form  to  the  above  and  can  be  written  down  by  inspection  using 
analogous  terminology.  To  wit,  when  =  1,  then: 


0  d9  /  I  \  .  j,  • 

"  J  dz*  H 

HO  “ 


(6.53) 


where  o  =  o,  ,/3.  If  a  <  o„,  then  the  response  is  elastic  and 
kk  -  I 


0  -  K  • 

o  kk 


(6.54) 


Of  course,  the  above  discussion  applies  only  to  the  case  where  is  equal 


to  unity. 
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7.  Elastic  Response  as  Points  of  Interdeterminancy  of 


dZ^ 


In  this  Section  we  shall  show  that  when  =  1,  there  exist  physical 
processes  for  which  is  equal  to  zero  and  that  in  fact  these  processes 
are  associated  with  elastic  deformation.  Analogous  conclusions  can  be  drawn 
with  regard  to  d?^  when  =  1.  The  above  remarks  are  made  in  the  context 
of  assumed  strict  independence  between  deviatoric  and  hydrostatic  response, 
in  the  sense  that  a  history  of  deviatoric  strain  has  no  effect  on  the  hydro¬ 
static  response  and  correspondingly  a  history  of  hydrostatic  strain  has  no 
effect  on  the  deviatoric  response.* 

It  is  important  for  our  purposes  to  introduce  certain  definitions: 
Deviatoric  plastic  strain  space:  The  space  of  deviatoric  strain  com¬ 
ponents  with  metric  6^^.  The  coordinates  of  this  space,  will  be  denoted  by 
X.,  (i  =  1,  2,  9)  corresponding  to  the  components  of  the  tensor  . 

Deviatoric  plastic  strain  path:  A  continuous  line  in  space. 
Remark:  The  path  determines  the  variation  of  x^  in  this  space.  We  limit 

the  discussion  to  paths  that  pass  through  the  origin.  In  this  event  =  0 
at  some  point  in  the  path,  in  fact  the  origin. 

Let  the  extremity  of  the  vector  x  change  continuously  along  the  path. 
We  define  quantities  s  and  ?  as  follows: 

s:  The  distance  of  the  point  x^  "the  origin  measiared  along  the 

*Henceforth  the  suffix  D  will  be  omitted.  Also,  this  chapter  will  apply 
exclusively  to  deviatoric  response  and  occasional  omission  of  the  word 
"deviatoric"  must  not  be  construed  to  imply  otherwise . 
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path,  where  X.  is  the  current  position  of  the  extremity  of  the  vector  x- 

C:  The  cummulative  distance  travelled  by  the  extremity  of  x»  i-e. 

d^  =  1 ds I  (7.1) 

or 

C  =  /|ds|  (7.2) 

Stipulations ; 

There  exists  a  point  on  the  path  at  which  s  =  0;  in  the  present  case 
this  is  the  origin  of  the  deviatoric  strain  space. 

There  exists  a  positive  direction  of  s ,  which  is  the  direction  of 
increasing  s  along  the  path.  In  a  similar  sense  there  exists  a  negative 
direction  of  s  which  is  the  direction  of  decreasing  s  along  the  path. 

The  path  is  separated  into  two  regions  of  negative  and  positive  s 
with  a  common  point  which  is  the  origin  of  x^»  s  is  equal  to  zero  at  this 
point . 

In  a  physical  context  we  note  that  in  the  case  of  simple  shear,  for 
instance,  the  shear  strain  is  positive  when  the  material  is  sheared  in  a 
specific  (positive)  direction  and  negative  when  shearing  takes  place  in 
the  opposite  direction.  In  a  similar  sense  in  a  uniaxial  plastic  strain 
experiment  there  exists  a  positive  (tensile)  strain  and  a  negative  (com¬ 
pressive)  strain.  Of  course  this  is  precisely  the  meaning  of  s  in  one 
dimensional  tests. 

Remark  I.  Givaia  path,  the  derivative  ^  experiences  a  jump  equal 
to  2  when  one  traverses  the  origin  in  the  positive  s  direction.  The  jump 
is  equal  to  -  2  when  the  origin  is  traversed  in  the  negative  s  direction. 
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The  proof  is  trivial:  3—  =  1  at  s  =  0+;  3—  =  -1  at  s  =  0-. 

d?  dc 

ds 

Remark  II.  The  derivative  -r—  is  indeterminate  at  s  =  Q  =  0. 
-  dc 

Figure  7  illustrates  the  above  statements. 

Discussion  of  eg.  (6.36)  at  ^  =  0. 


At  c  =  0,  and  therefore  at  s  =  0 


s 


s 


Y  dz 


Therefore , 


since  at  ?  =  0,  f(o)  =  1. 

On  the  basis  of  eq.  (7.4) 


(7.3) 


(7.4) 


(s 


^.2 

Ydc^ 


(7.5) 


But  the  right  hand  side  of  eq.  (7.5)  is  indeterminate  at  ?  =  0,  hence 
||s||,  and  therefore  s,  are  also  indeterminate.  However  C=  0,  requires 
that 

=  0  (7.6) 

Equation  (7.6),  therefore,  and  the  definition  of  (eq.  (6.5)  with  =  1) 
lead  to  the  relation, 

s  =  2vie  (7.7) 

which  is  a  statement  to  the  fact  that  at  ^  =  0,  the  deviatoric  response  is 


elastic . 
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Eqviation  (6 .36)  at  C  =  0+  0-. 

In  the  vicinity  of  s  =  0,  let ^  be  given  by  the  relation 

=  £s  (7.8) 

where  .  are  the  direction  cosines  at  s  =  0.  Then  as  a  result  of  eq.  (6.36) 
ID 

(note  that  the  second  term  of  the  right  hand  side  of  eq.  (6.36)  tends  to 
zero  at  z  ->  0 ) 

s  .  =  Sv^-.  .  (7.9a) 

s  =  -s°Jl..  (7.9b) 

where  s^  =  s(0+)  and  s  =  s(O-).  The  de viator ic  stress  tensor  Is  therefore 
discontinuous,  with  a  discontinuity  of  magnitude  2sy,  at  s  =  C  =0. 

Equation  (6.36)  at  Reversal  Points. 

Definition:  A  reversal  point  R  is  a  point  on  the  strain  path  s  at 
which  there  is  a  reversal  in  the  sign  of  the  strain  increment.  In  effect  at 
R 

dx.  =  -  dx."*"  (7.10) 

1  1 

With  reference  to  Figure  8,  the  following  relations  are  of  interest: 

dx- 

(R')  =  (R)  +  0.  (7.11) 

ds  ds  1 

Lim  9.  =  0,  R'  R. 

1 
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However,  as  a  result  of  our  definitions  of  ds  and  dC: 

Li.  I  I,,  =  1 
R'  R 


Lim 

R""  R 


Hence,  invoking  eq.  (6.10): 


I'i’"  ^  1  . 

R^  ^  R  dz  R 


fCOiR 


Lim 

R"  R 


ds 

dz 


-  «^>Ir 


(7.12a) 


(7.12b) 


(7.13a) 


(7.13b) 


Also, 

Lim  r(R')  =  r(R)  (7.14a) 

R"  -5"  R 


Lim  r(R"')  =  r(R) 

R""  ^  R. 

Let  in  the  vicinity  of  R 

d/^  =  Ms  (7.15) 

where  are  the  direction  cosines  of  the  tangent  to  the  path  at  R.  Then 
as  a  result  of  eq. ’s  (7.11)  throu^  (7.15)  and  eq.  (6.36): 
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s(R')  =  +  r(R),  (7.16a) 

s(R"')  =  -  s°Jl.  .fp(^)  +  r(R),  (7.16b) 

in  the  limits  R^  ->■  R,  R"  R. 

Since  in  fact  is  normal  to  the  stress  hypersphere  at  R  and 
SYf(C)  is  the  current  radius  of  the  hypersphere,  eq.'s  (7.16a,b)  admit  the 
geometric  construction  shown  in  Figure  9 . 

Note  that  the  response  in  the  stress  range  R''R''‘'  is  elastic,  and 
corresponds  to  constant  ?  and  z.  During  the  elastic  response  d?  is  there¬ 
fore  zero  and  the  constitutive  equation  is  simply 

ds  =  2]j  de  (7.17) 

~  o  , 

A  plot  of  C  versus  s  for  this  process  is  shown  in  Figiire  10. 

ds  • 

Clearly,  at  R,  the  slope  ^  is  indeterminate. 


More  General  Deformation  Processes  With  Constant 


We  begin  with  eq.'s  (6.36)  and  (6.37)  which  we  write  in  the  form 

-o  ~  +  r  (7.18) 


s  =  s 


Y  dz 


or,  equivalently j 


s  = 


s?/ 

0 


-  S'  "  'I 


(7.19) 
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Sinc6,  however,  the  delta  function  in  eq.  (7.19)  was  obtained  by 
a  limiting  process  in  which  ->  1,  we  replace  6(  z  )  by  a  fxmction  which 
in  the  above  limiting  process  becomes  a  delta  function,  i.e.,  we  write 

2u  z  ,  d^ 

o  -  —2.  f  2  ^  ^  dz"  +  r(  z  )  (7.20) 

s  -  T__v_  J  «  dz'  . 


where 

2y 

=  o  .  (7.21) 

(l-kj^)s° 

and  s  is  obtained  from  eq.  (7.20)  by  a  limiting  process  in  which  tends 
to  unity.  Equation  (7.20)  may  now  be  written  as  a  differential  equation 
in  the  form : 


1-k 

—  (s  -  r)dz  +  Tj — —  <i(s  -  r)  =  d^  . 


(7.22) 


We  note  in  passing  that  in  the  limit  k^  1 


s  -  r 


=  S:jf(?)  =  Sy 


(7.23a,  b) 


a  result  given  previously  in  eq.  (6.41), 

d?  : 


We  also  recall  the  definition  of 


dZ  =  de  -  ds 


(7.24) 


Equations  (7.22)  and  (7.24)  combine  to  give  the  result 


k-  —  dc  +  (1-k  )de  =  dZ 
Is  i  ~  ~ 

Y 


(7.25) 
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where 


s  =  s  -  r 


de  =  de  -  dr 


(7.26) 
(7  J27) 


However,  from  eq  (6.37) 


dr  =  2p  (p, (o)dU  +  h  dz) 
^  o  1  ?  ~ 


where 


h  =  f  P'( 


d? 


(7.28) 


(7.29) 


and 


Pl  =  -dl  (0  <  z  <.) 


(7.30) 


Thus,  p"  is  well  behaved  (continuous,  monotonically  decreasing)  function  of 
z.  Note  that  h(0)  =  0. 

Equations  (7.25)  and  (7.28)  combine  to  give  the  following  equation 

for  dz: 

adz  +  (l-k^)de  =  cd^ 


(7.31) 


where 


a  =  {k,  ^  -  k.(l  -  k  )h}/f(c) 
1  O  1  1  ^ 

Sy 


c  =  1  +  k^(l  -  k^)p^(o) 


(7.32) 


(7.33) 


By  taking  the  norm  of  both  sides  of  eq.  (7.31)  one  obtains  the  following 
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quadratic  equation  in  d  ? 


2  a  •  ded?  -  (l-k^)l lde| =  0  (7.34) 


The  two  roots  of  eq.  (7.34)  are  given  by  eq.  (7.35): 


d? 


a  •  de  ±  /  a  •  de  +  (l-k^)C| |de|  P 
C 


(7.35) 


where 


(7.36) 


It  may  be  shown  that 


Lim  C  = 


k^-.! 


where 


In  the  limit  of  k^  ->  1,  the  two  roots  of  dC  are,  therefore: 


(7.37) 


(7.38) 


(7.39) 


with  the  constraint  d  C  ^  0 . 
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The  Two  Roots  ofdC- 


At  the  initiation  of  plastic  deformation  h 
denced  by  eq.  (7.29).  Now,  from  eq.’s  (6.29)  and 


is  equal  to  zero  as  evi- 
(6.30) 


n-1 

I 


r=l 


R  e  ^r 
r 


(7.40) 


where  R  are  all  positive  as  shown  in  Appendix  B*  Therefore,  p^(o)  is  posi- 
tive.  However,  the  scalar  product  s*h  may  be  negative.  In  fact,  in  the 
particular  case  where  consists  of  a  single  exponential  term,  i.e. 

P^(C)  =  fR^e'^l^  (7.41) 


The  function  p'(c)  is  given  by  eq.  (7.42): 


p-(C)  = 


-3  C 


(7.42) 


in  which  event 


h(c)  =  -6 


r(?) 

1  w 


Thus 


s  ‘h 


s  ‘r 


(7.43) 


The  scalar  product  s*h  will,  therefore,  be  negative  if  the  angle  between 

A 

s  and  r  is  acute  —in  this  particular  case.  It  is  therefore  quite  possible 
that  during  the  process  of  deformation  will  decrease  from  a  positive 


value  to  zero. 
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We  shall  first  investigate  the  situation  where  is  bounded  by 
the  inequalities 

0  <  <  2{1  +  p^(0)} 

Case  (i):  s*de  >  0 

In  this  case  the  dissipation  is  positive,  see  eq.  (2.14),  since  dC 
is  in  fact  positive  and  is  given  by  the  eq.  (7.44) 

A 

s*de 


d? 


(7.44) 


The  process  is  admissible. 


Case  (ii):  s*de  <  0 

In  this  case  the  dissipation  is  negative  (since  dc  is  negative, 
see  eq.  (2.14))  and  the  process  is  inadmissible.  Therefore  the  second 
root  of  d^  must  be  chosen  and  as  a  result : 

dC  =  0.  (7.45) 

In  this  event  the  dissipation  is  zero  and  the  deformation  is  elastic 
since  eq.  (7.45)  implies  that 

d<^=  0  (7-46) 

in  which  case 

ds  =  2u  de.  (7.47) 

o  ~ 

Remark;  Case  (ii)  is  the  "unloading  condition"  given  unequivocably  by  the 
endochronic  theory.  More  specifically,  given  a  state  of  stress  on  the 


V 
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y£^02.d.  suipfsic©^  irli0  pi?oc0ss  is  ciissip3.triv©  3.tic1  ^  Oj  if  "the  S'ti'Sin  incr© 

ment  makes  an  acute  angle  with  the  radius  vector  through  the  point  which 
represents  the  state  of  stress;  the  process  is  elastic  if  the  angle  is  equal 
to  or  greater  than  ninety  degrees. 

The  Case:  =  0 

The  situation  where  ->  0  results  in  d?  becoming  unboxmdedly  large 
in  which  event  the  material  flows  without  limit .  This  event  will  occur 
when 

^h  =  -  s^d  +  Pi(o))f(i;)^  (7.48) 

Stress  Paths  Within  the  Yield  Surface. 

If  eq.  (6.36)  is  to  apply,  dC  must  be  different  from  zero  (in  fact 

d'2 

positive),  since  otherwise  the  derivative  not  determinate.  There™ 

fore,  if  eq.  (6.36)  applies,  then  the  incremental  process  emanating  from 
state  s  is  dissipative  and  vice-versa.  In  the  same  vein,  eq.  (6.51)  is  a 
direct  consequence  of  eq.  (6.36).  Therefore  if  an  incremental  process  ori¬ 
ginating  at  state  s  is  to  be  dissipative  eq.  (6.51)  must  necessarily  apply. 
It  follows  as  a  corollary  that  stress  states  which  do  not  satisfy  eq.  (6.51) 
cannot  lead  immediately  to  dissipative  processes. 

Remark ;  State  s  which  lie  within  the  yield  surface  cannot  lead  immediately 
to  dissipative  Incremental  processes. 

It  follows  from  the  above  remark  that  all  incremental  processes 
which  emanate  from  stress  states  within  the  yield  surface  are  elastic 
(dC  =  0),  in  which  event  for  all  ds  emanating  from  such  state  s: 


ds  =  2yde 


(7.49) 
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Incremental  Determination  of  the  Constitutive  Response. 

Given  that  the  constitutive  properties  of  a  material  are  known, 
let,  at  some  point  in  the  course  of  the  deformation  (strain)  history, 
s  and  h  also  be  known.  Then  can  be  determined  from  eq.  (7.38). 

Given  an  increment  of  strain  de,  eq.  (7.39)  determines  d^,  depending  on 
the  sign  of  s*de.  If  the  sign  is  negative  or  zero,  then  d?  =  0  and  ds 
is  given  by  eq.  (7.47).  If  the  sign  is  positive  then 


,  _  d? 

■  f(0 


(7.51) 


In  the  limit  of  k^^  1,  then  as  a  result  of  eq.  (7.22)  or  eq.  (7.18) 


A 

dl?  =  —  dz 

I*  o 

Hence,  di?  may  be  calculated.  Knowledge  of  dj  leads  to  the  direct  deter¬ 
mination  of  ds  from  eq.  (7.24)  in  the  limit  of  k^  ^  1,  i.e., 

ds  =  2v  (de  -  d?) 

o  - 

Clearly  dr  may  now  be  found  from  eq.  (7.28),  and  the  procedure  may  thus  be 
repeated  at  will. 

We  note  that  at  the  initiation  of  plastic  deformation  ^  ~  ^  ^  2 

and  s  =  2y  e. 

o_ 
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8.  The  Closure  of  Hysteresis  Loops. 

The  question  of  closure  is  treated  by  means  of  an  application  to 
aluminum  undergoing  stress  reversals  in  a  uniaxial  stress  field.  One  has 
difficulties  in  obtaining  closed  hysteresis  loops  in  the  first  quadrant  of 
the  stress-strain  diagram,  if  one  uses  the  definition  of  intrinsic  time  of 
Reference  1,  i.e.,  eq.  (1.1).  The  purpose  of  this  Section  is  to  demonstrate 
that  no  such  difficulties  arise  when 


dC  =  |d0j 


and 


(8.1) 


(8.2) 


and  k  =  1. 

For  the  purpose  of  demonstration  we  use  the  constitutive  equation 


(8.3) 


where 

E(C)  =  +  Ej 

2 

with  the  following  values  of  the  constants:  =  5  tons/in  ,  E^  =  1500 

tons/in^,  E^  =  200  tons/in^,  a  =  500.  The  elastic  modulus  E^  is  equal  to 
4.46  X  10^  tons/in^l  however,  this  value  is  not  pertinent  to  the  present 
application  since  is  plotted  directly  against  0^,  The  computation  of 
the  stress  response  to  loading-unloading-reloading  histories  is  straight¬ 
forward  and  will  not  be  treated  in  detail  here.  For  a  more  elaborate  treat 
ment  of  this  problem  the  reader  is  referred  to  Reference  2,  The  results 
are  shown  in  Figure  ii. 
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The  point  A  is  the  initial  yield  event.  Point  B  denotes  the  point  of 
unloading  and  C,  D,  E,  F,  G,  and  H  are  points  of  reloading.  Note  that  points 
such  as  C,  where  o  >  -  20^,  gives  rise  to  reloading  paths  which  coincide 

with  the  unloading  paths.  In  such  a  case  the  area  of  the  hysteresis  loops  is 
zero.  However,  reloading  from  points  where  a  <  -  2a^  gives  rise  to  loops 

that  are  closed,  as  shown  in  the  Figure.  In  this  constitutive  theory  as  it 
is  constituted  by  equations  (7.1),  (7.2)  and  (7.3),  no  open  loops  exist  for 
k  =  1. 

9.  Discussion 

In  this  paper  we  have  presented  an  endochronic  theory  of  plasticity 
which  is  based  on  a  new  definition  of  intrinsic  time  which  is  broader  than 
the  one  we  proposed  before. In  its  two  limiting  forms  the  intrinsic  time 
measure  is  expressible  either  in  terms  of  the  increment  of  the  total  defor¬ 
mation  on  one  hand,  or  in  terms  of  the  increment  of  the  plastic  deformation, 
on  the  other.  An  intermediate  definition  is  made  possible  by  means  of  a 
material  tensor 

The  theory  is  shown  to  predict  essential  featvires  of  the  plastic  re¬ 
sponse  of  metals  and  in  certain  asymptotic  cases  it  evinces  commonalities 
with  classical  theories  of  plasticity  and  contains  a  number  of  classical  plas- 
"ticity  models  as  special  cases.  In  particular,  the  existence  of  a  yield  sur¬ 
face  is  shown  as  a  consequence,  when  the  intrinsic  time  measure  is  defined  in 
terms  of  the  increment  of  plastic  strain.  Also,  the  response  of  an  elastic 
perfectly  plastic  solid  is  shown  to  be  predictable  by  the  theory.  However, 
the  most  significant  feature  of  the  present  theory  is  that  it  predicts  correctly 
in  the  case  of  metals  where  is  almost  equal  to  unity  the  stress  response 
to  strain  histories  involving  reversals  in  the  sign  of  the  strain  rate.  This 
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is  a  particularly  powerful  result,  as  one  and  the  same  constitutive  equation 
is  used  in  the  prediction  of  the  stress  response  to  such  reversals.  When  k 
is  equal  to  unity  a  dichotomy  in  the  constitutive  representation  results  in 
the  sense  that  an  elastic  constitutive  relation  must  be  used  in  the  course  of 
strain  histories  involving  "unloading" . 

It  is  felt  that  the  endochronic  theory  has  now  evolved  to  a  stage 
where  it  can  be  used  in  design  as  a  predictive  tool  in  the  understanding  of 
the  mechanical  response  of  metals  and  other  materials. 
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Appendix  ^ 


We  solve  eq.  (6.19)  by  the  Laplace  Transform  technique,  where 


-a  2  -a  2 

6(2)  =  e  +626 


(A.l) 


Then,  as  a  result  of  eq.  (6,19), 


P  = 


p  4-  a 


(1-k^)  p  +  (o^+  Oj-  ka)  p  + 


(A. 2) 


where 


=  ®2“i  +  <^“2  * 


(A.3) 


It  can  be  shown  by  straightforward  analysis  that  as  k^^  1 


.£_L_a 


n-k^Kp+e^TTp+i^ 


(A. 4) 


where 


6.  = 


(a,G.  +  a-G.) 


1  1-k^  'll  22 


(A. 5) 


So  = 


“l“2 


2 


(A. 6) 


In  the  limit  of  k^-»-  1  and  as  a  result  of  eq.  (A. 4),  p(2)  may  be  written 
explicitly  in  the  form; 
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where 

It  is 
of  eq 


“l“2 


(A. 7) 


X  =  Gj_a^  +  . 


(A. 8) 


worth  repeating  that  as  1,  the  first  term  on  the  right-hand  side 

,  (A. 7)  becomes  proportional  to  a  delta  function  and 

/o«-  a,  \2  o.a. 

p(z)  =  Y  >5(2)  +  G-,G-  (  — -)  expi--^z).  (A. 9) 


The  fact  that  6(z)  is  an  essential  part  of  p(z)  when  k^^  =1,  was  demonstrated 
in  the  text  in  the  case  of  any  number  of  internal  variables . 
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Appendix  ^ 


We  note  the  following  relations: 


R  =  Q  -  pP 


(B  -1) 


P  =  Q  I 


(B-2) 


P  +  “i 


We  also  note  that,  since  -3^  are 


the  zeros  of  R,  it  follows  from  eq.  (B-1) 


(B-3) 


Also  from  eq.  (B-2): 


«-Bs>  =  o'-Ss*  I 

r  s  r 


(B-4) 


It  follows  from  eq. ’s  (B-3)  and  (B-4)  that 


r  r  _  _  _±. 

^  -3  +  a  3„ 

r  s  r  s 


(B-5) 


Also  from  eq. ’s  (B-1)  and  (B-2) 


pG 

1  -  I  - ^ 

r  ^  “r  p  =  -3, 


CB-6) 


But,  from  eq.  (6.22), 


I  G  =  1 

L  r" 


(B-7) 


Hence  from  eq.'s  (B-6)  and  (B-7) 
ot  G 

y  —  =  1 

^  -3  +  a 

r  s  r 


(B-8) 
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At  this  point  we  are  in  a  position  to  evaluate  R^.  To  this  end  we  note  that 
as  a  result  of  eq.  (6.24), 


=  Q  y  — 

^  ^  ^  D  +  a 


(B-9) 


Also  in  view  of  eq.  (B-2): 


=  Q I  TTir  iiTfV  - ^  (p  A  )2 


Cb-10) 


r  r 


Therefore  It  follows  from  eq.  (B-1)  and  eq.’s  (B-9)  and  (B-10)  that. 


Iftit 

T  ^  r  r  r  Ir  rr  r 


^r  ) 

■  TpT^)2j 


(B-ll) 


At  this  point  evaluating  at  -B  and  invoking  eq.  (B-6),  we  find  that 

a 


(-3s)  =  -Q(-3s)  I  +  I 

■  ^  g  ^  r  s  r 


(B-12) 


r(.6  )  =  -Q(-e^)  I  \—y, 
r  r  s 


(B«13) 


Now  from  (B-1) 


Q(-3s) 
P(-33)  =  ^ 


CB-14) 


At  this  point  use  of  eq.  (6.28)  gives  an  explicit  expression  for  given  by 
eq.  (B-15) 


=  6 


(B-15) 


49 


However,  the  constants  a^,  3^  and  are  all  positive.  Thus  the  residues 
are  all  positive.' 
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Appendix  C 

In  the  theory  of  Laplace  Transform,  one  may  deal  with  fvinctions  with 
simple  jump  discontinuities  by  enlarging  the  class  of  differentiable  functions 
to  include  such  functions.  This  is  achieved  by  making  the  Heavinside  step 
function  a  member  of  this  class  and  by  stipulating  that  the  derivative  of  this 
function  is  the  Dirac  delta  function,  i.e., 

|^H(t)  =  6(t)  (C-1) 

If  one  does  not  do  this  one  can  only  deal  with  differentiable  functions  in  the 
classical  sense  in  which  case  if  such  a  function  f(t)  is  such  that  f(0)  ^  0, 
then  f(t)  is  defined  in  the  interval  0  <  t  <  «  and,  of  course,  as  is  well 
known 

l(^)=  pf(p)  -  f(0+).  (C-2) 

where  L  denotes  the  Laplace  transform  over  the  open  interval  (0, «  )  with  respect 

to  the  parameter  p. 

It  follows  from  (C-2)  that 

pf(p)  =  f(0+)  +  L  (^)  (C-3) 

or 

L‘V(p)  =  f(0+)6(t)  +  ^  I  (C-4) 

If,  however,  functions  with  simple  jump  discontinuities  are  admitted,  the 
Laplace  transformation  may  be  used  over  the  half  open  interval  )  and  a  func¬ 
tion  can  now  be  defined  on  this  interval  with  the  condition  that  f(0)  =  0, 

f  (t)  =  f(t),  t>0.  This  is  what  we  have  done  in  this  paper.  This  function 
o 

f  (t)  now  has  a  jump  discontinuity  at  t  =  0  of  strength  f(0+). 
o 
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In  this  case 

Note  that  f  (p)  =  f(p)  since  differs  from  f  only  on  a  set  of  measxare  zero. 
Thus,  using  eq.  (C-5) 

L  — 2-  =  pf(p)  .  (C-6) 

dt 

since  f  (o)  =  0. 
o 

It  follows  that 

df  jjT 

L”^pf(p)  =  =  f(0+)6  (t)  +  ^  1^^^  (C-7) 

in  accordance  with  eq.  (C-1). 
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Appendix  D 


Experimental  determination  of  the  functions  Pj^(z)  and 


In  pure  shear  eq.  (6.36)  becomes 


T 


o  dY*" 
’Y  dz 


+  2u 


z 

/p, 


(z  - 


(D-1) 


where  x  is  the  shear  stress  and  Y^  the  plastic  shear  strain, 
relation 


dz 


- 

"  f(0 


We  recall  the 

(D-2) 


The  suffix  D  on  C  and  z  has  been  omitted  for  simplicity  of  notation . 

With  reference  to  Figure  12  which  is  self-explanatory  and  eq.  (D-1), 
the  following  resiUts  are  evident: 

(i)  5  is  the  plastic  shear  strain  at  point  A,  and  is  equal  to  the 

distance  oA' . 

(ii)  The  equation  below  holds 


=■ 


(D-3) 


Thus  f ( 0  can  be  determined  from  measurements  of  the  yield  points  upon  load¬ 
ing  (point  Y)  and  load  reversal  (point  Y').  The  function  z(0  may  now  be 
determined  by  integration  of  eq.  (D-2).  Thus  we  have  the  functional  relation 


f(0  =  f{?(z)>  =  f(z) 

At  this  jimcture  we  show  that  p^(z)  can  be  determined  from  a  simple  monotonic 

loading  test.  During  such  a  test  dz  is  related  to  dy  by  eq.  (D-5) 

.  - 
"  fCyP) 


(D-5) 
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since,  in  fact,  during  loading 


d?  =  dy^  (D-6) 

and 

5  =  yP  (D-7) 

upon  observation  of  the  condition  ^  =  0  whenY=  0.  As  a  result  of  eq.'s 
(D-4)  through  (D-7),  eq.  (D-1)  reduces  to  the  following  form: 

Z 

T  =  s^fCz)  +  2y  /  p^(z  -  z'')f(z‘')dz''  (^'8) 

^  °  o 

However  since  t(y^)  is  known,  then  so  is  and  hence  tCz)  =  tCc(z)}. 

Thus,  since  f(z)  is  also  known,  the  function  g(z),  where 


g(z)  =  t(z)  -  SyfCz) 

is  a  known  function  of  z.  Equation  (D-8)  therefore,  now  reduces  to  the  fol 
lowing  Volterra  integral  equation 

2v  /  f(z  -  z')  p  (z')dz'  =  g(z) 


wh.6r*s  "tlis  nghti  li3.Tid  sid.©  is  3.  known  function# 

Equation  (D-10)  determines  p^(z)  uniquely  if  f  is  a  continuous  bounded 
function  of  z  in  the  interval  (0,z).  The  function  f  is  expected  to  have 
these  properties  on  physical  grounds. 


References 


54 


1.  Valanis,  K.  C.,  "A  Theory  of  Viscoplasticity  without  a  Yield  Sxjrface, 

Part  I:  General  Theory,"  Archives  of  Mechanics,  Vol.  23,  p.  517  (1971). 

2.  Valanis,  K.  C.,  "A  Theory  of  Viscoplasticity  without  a  Yield  Surface, 

Part  II:  Application  to  Mechanical  Behavior  of  Metals,"  Archives  of 
Mechanics ,  Vol.  23,  p.  535  (1971). 

3.  Valanis,  K.  C.,  "Effect  of  Prior  Deformation  on  the  Cyclic  Response  of 

Metals,"  J .  App .  Mech . ,  Vol.  41,  p.  441  (1974). 

4.  Valanis,  K.  C.  and  Wu,  H.  C. ,  "Endochronic  Representation  of  Cyclic  Creep 

and  Relaxation  of  Metals,"  J.  App.  Mech.,  Vol.  42,  p.  67  (1975). 

5.  Valanis,  K.  C.,  "On  the  Foundations  of  the  Endochronic  Theory  of  Plasti¬ 
city,"  Archives  of  Mechanics,  Vol.  27,  p.  857  (1975). 

6.  Bazant,  Z.  P.,  "A  New  Approach  to  Inelasticity  and  Failure  of  Concrete 
Sand  and  Rock:  Endochronic  Theory,"  Proceedings,  Society  of  Engineering 
Science,  11th  Annual  Meeting,  Durham,  North  Carolina  (1974). 

7.  Bazant,  A.  P.,  "Endochronic  Theory  of  Inelasticity  and  Failure  of  Con¬ 
crete,"  Journal  of  Engineering  Mechanics  Division,  ASCE,  Vol.  102,  p. 
12360  (1976). 

8.  Wempner,  G.  and  Aberson,  J.,  "A  Formulation  of  Inelasticity  from  Thermal 
and  Mechanical  Concepts,"  Int.  Journal  of  Solids  and  Structures,  Vol. 

12,  p.  705  (1976). 

9.  Miller,  A.,  "Inelastic  Constitutive  Model  for  Monotonic,  Cyclic,  and 
Creep  Deformation.  1.  Equations  Development  and  Analytic  Procedures," 

J.  Eng.  Mater  Technol  Trans,  ASME,  Vol.  98,  No.  2,  pp.  97-105  (1976). 

10.  Hart,  E.  W. ,  "Constitutive  Relations  for  the  Nonelastic  Deformation  of 
Metals,"  J.  Eng.  Mater  Technol  Trans,  ASME,  Vol.  98,  No.  3,  pp.  193- 
202  (1976). 

11.  Bodner,  S.  R. ,  "Constitutive  Equations  for  Elastic-Viscoplastic  Strain- 
Hardening  Materials,"  J .  App .  Mech . ,  Vol.  42,  No.  2,  pp.  385-389  (1975). 

12.  Havner,  Kerry  S.,  "Some  Concepts  and  Results  from  the  Mechanics  of 
Crystalline  Solids  ...,"  Proceedings,  Workshop  on  Applied  Thermovisco- 
plasticity.  Northwestern  University  (October  1975). 

13.  Valanis,  K.  C.,  "Irreversibility  and  Existence  of  Entropy,"  Journal 
of  Non-Lin.  Mech.,  Vol.  6,  p.  337  (1971). 

14.  Krempl,  E.,  "Plastic  Deformation,  Microstructure  Change  and  the  Variable 
Heredity  Property  of  Metals,"  Proceedings,  Workshop  on  Inelastic  Consti- 
tutive  Equations  for  Metals,  Rensselaer  Polytechnic  Institute,  Troy, 

New  York  (1974). 


59 


Fig,  5,  Dissipation  in  loading  and  unloading, 


Fig.  11.  Theoretically  predicted  unloading -reloading  behavior  of  aluminum. 


